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Disc-jet coupling in the 2009 outburst of the black hole 
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ABSTRACT 

Wc present an intensive radio and X-ray monitoring campaign on the 2009 outburst 
of the Galactic black hole candidate X-ray binary H1743-322. With the high angular 
resolution of the Very Long Baseline Array, we resolve the jet ejection event and mea- 
sure the proper motions of the jet ejecta relative to the position of the compact core 
jets detected at the beginning of the outburst. This allows us to accurately couple the 
moment when the jet ejection event occurred with X-ray spectral and timing signa- 
k><| \ tures. We find that X-ray timing signatures are the best diagnostic of the jet ejection 

. event in this outburst, which occurred as the X-ray variability began to decrease and 

' the Type C quasi-periodic oscillations disappeared from the X-ray power density spec- 

trum. However, this sequence of events does not appear to be replicated in all black 
hole X-ray binary outbursts, even within an individual source. In our observations of 
H1743-322, the ejection was contemporaneous with a quenching of the radio emission, 
prior to the start of the major radio flare. This contradicts previous assumptions that 
the onset of the radio flare marks the moment of ejection. The jet speed appears to 
vary between outbursts, with a possible positive correlation with outburst luminosity. 
The compact core radio jet reactivated on transition to the hard intermediate state at 
the end of the outburst, and not when the source reached the low hard spectral state. 
Comparison with the known near-infrared behaviour of the compact jets suggests a 
gradual evolution of the compact jet power over a few days near the beginning and 
end of an outburst. 

Key words: accretion, accretion discs - black hole physics - ISM: jets and outflows 
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1 INTRODUCTION 

The nature of the coupling between inflow and outflow in 
accreting compact objects has been the focus of intensive 
study in recent years. Galactic X-ray binary sources pro- 
vide ideal laboratories in which to probe this disc-jet cou- 
pling, since they pass through a number of well-defined ac- 
cretion regimes as they evolve through their full duty cycles 
on timescales of months to years. Thus we can study a range 
of distinct accretion states in the same object, hence control- 
ling for the effects of compact object mass, donor type and 
environment on the formation of jets in different accretion 
regimes. 

The properties of the accretion flow are best studied 
in the X-ray band, and have b een categorise d into a num- 
ber of canonical states (see, e.g. lBellonill2010l . for a review), 
which are identified by specific X-ray spectral and timing 
characteristics. Black hole X-ray binaries spend the major- 
ity of their time in a very low luminosity (< 1 Q 335 ergs -1 ) 
quiescent state (|McClintock fc Remillardll2006l '). At the be- 
ginning and end of an outburst, such systems are observed 
in the low/hard state (LHS), with an X-ray spectrum dom- 
inated by a fairly hard power law component (photon index 
r ~ 1.5), and a power density spectrum showing strong 
band-limited noise with integrated fractional rms variability 
~ 0.3. During the peak phase of an outburst, a black hole 
X-ray binary will make a transition through a set of inter- 
mediate states to the high/soft state (HSS). The HSS has 
an X-ray spectrum dominated by soft thermal blackbody 
emission from a geometrically thin, optically thick accre- 
tion disc, together with a weak, steep (r > 2.5) power- 
law component, and showing a power density spectrum 
with a much reduced fractional rms va riability of < 0.1. 
The in termediate states were classified bv lHoman fc Bellonil 
(|2005h into two distinct classes, the hard and soft interme- 
diate states (HIMS and SIMS respectively), distinguished 
by both the photon index of the power law component and 
by their X-ray variability properties, notably a rapid de- 
crease in the fractional rms variability and a transition be- 
tween two types of quasi-periodic oscillations (QPOs) on 
moving from the HIMS to the SIMS. At the end of the 
outburst, the X-ray luminosity decreases, and the source 
makes a reverse transition from the HSS through the HIMS 
back to the LHS. The reduced luminosity of the soft-to-hard 
transition as compared to the hard-to-soft transition leads 
to an observed hysteresis in a hardness-intensity diagram 
(HID). While there ar e alternative sets of state definit ions, 
most notably those of iMcClintock fe Remillardl (|2006l ). we 
have chosen to use the terminology outlined above pri- 
marily for consistency with much of the published liter- 
ature on the coupling of radio and X-ray emission from 
black hole X-ray binar ies, specifically the phenomeno logi- 
cal picture presented bylFender. Belloni fe Galld (|2004T ) and 
IFender, Homan fc Bellonil (|2009l V 

The relativistic jets in black hole X-ray binaries are 
seen in two distinct forms. Steady, compact jets have 
been directly imaged in hard X-ray spectral states (e.g . 
iDhawan. Mirabel fc Rodriguez! l200ol ; IStirling et al] l200lf ). 
and are o bserved to h ave a flat or slightly inverted radio 
spectrum (iFendeitoOll ). They are best modelled as partially 
self-a bsorbed conical outflows (e.g. Blandfor d fc Konigll 
1 19791 ). and are believed to carry away a significant, if not 



dominant , fraction of the power liberated by the accre- 
tion flow (|Fender. Gallo fc Jonker1l2003l ; iGallo et al.l 120051 ). 
The second, more spectacular, type of jet morphology is 
seen during outbursts, as bright, steep-spectrum, discrete, 
relativistically-moving knots of radio emission. In most sys- 
tems, these are only observed on milliarcsecon d scales using 
very l ong baseline interferometr y (VLBI; e.g. Tingav et al.l 
1 19951 ; iHiellming fc Rupenl 1 19951 ; iHannikainen et all |200]T k 



although for the brightest systems, they may be ob- 
served on larger angular scales with connected-element 
int erferometers such as the Very Large Array (VLA; 
e.g.|Mirabel fc Rodri guez 1994; M artf, Paredes fc Peracaulal 
l200lK ~ Optically thin radio emission may also be observed 
some distance downstream from the X-ray binary itself if 
the jets run into a particularly dense patch of the surround- 
ing interstellar medium, causing particle acceleration at the 
sites of the resulting external shocks (e.g. lCorbel et al"]|2002l . 
2005). Since this is dependent on the properties of the sur- 
rounding medium rather than the launching mechanism of 
the jets themselves, we will not focus on it in any detail. 

The relationship between the radio emission (assumed 
to arise from the jets) and the X-ray e mission from the ac - 
cretion flow was considered in detail bv lCorbel et~aH l|2004l) . 
who used the Australia Telescope Compact Array (ATCA) 
to sample the radio emission throughout the 2001 outburst 
of the black hole X-ray binary XTE J1650-500. Compar- 
ing their data to similar outbursts in three other systems, 
they found that the hard state (corresponding to the LHS 
in our terminology) is associated with emission from a par- 
tially self-absorbed compact jet, which persists into the in- 
termediate state (equivalent to the HIMS). Ejection events 
are associated with transitions from the intermediate to the 
steep power law state (HIMS to SIMS transition), and the 
jet is suppressed in the thermal dominant state (HSS), al- 
though residual radio emission may be observed as a result 
of interactions between the ejecta and the surrounding en- 
vironment. 

IFender fc B elloni ( 2 004h also made a detailed study of 
the relationship between the X-ray and radio emission (the 
disc-jet coupling) in a black hole X-ray binary, in this case 
focussing on the system GRS 1915+105. Although the X-ray 
spectral states in GRS 1915 +105 are non- s tanda rd, multi- 
wavelength observations by iMirabel et al. I l|l99St ) also pro- 
vided evidence that jet ejection events are triggered at the 
end of a phase of X-ray spectral softening. This was con- 
firmed by high-reso lution observations of a large radio flare 
jFender et al.ll 19991 ') . which measured the proper motions of 
the relativistic jets, showing them to be ejected as the source 
was moving from a relatively hard 'plateau' state to a softer 
X-ray state. 

IFender et alj (|2004T ) generalised this picture by synthe- 
sizing the existing data on the disc-jet coupling in several dif- 
ferent black hole X-ray binaries into a 'unified model'. In its 
simplest form, this model holds that the steady, compact jets 
exist during the quiescent state and LHS, with a relatively 
low bulk velocity (Lorentz factor T < 2). As the X-ray lumi- 
nosity (I/x) rises at the beginning of an outburst, so too does 
the jet power and hence the radio lumi nosity (L r ), follow- 
ing the non-linear correlation L r oc L^C (|Corbel et al.l [20o3 : 
IGallo. Fender fc Pooley||2003h . As the source moves into the 
intermediate states, the X-ray spectrum softens and the jet 
Lorentz factor increases, with the jet becoming unstable. At 



© 2012 RAS, MNRAS 000.1111201 



Disc-jet coupling in EVi %3-322 3 



a certain point in the HID, dubbed the 'jet-line', the core jet 
switches off, and a radio flare may occur as the faster mov- 
ing material ejected more recently catches up with the pre- 
existing, slower-moving flow, giving rise to internal shocks, 
which accelerate particles to produce the observed bright, 
optically thin, relativistically-moving ejecta. Alternatively, 
the shocks could be caused by increased outflow flux at con- 
stant velocity, rather than by a change in the bulk veloc- 
ity. Supporting such a scenario, the disappearance of the 
non-thermal power-law component of the X-ray spectrum 
at this time has led to suggestions that the shocks could 
be caused by the ejection of the corona and its subsequent 
interaction with material from the pre-existing steady jets 
jRodriguez. Corbel fc Tomsicld [20031 ; IVadawale et alj|2003t 
iRodriguez et al.l 120081 ). If the X-ray behaviour oscillates 
back and forth across the jet line, multiple radio flares may 
be observed, with the compact core jet reforming each time 
the spectrum becomes harder, and another ejection event oc- 
curring at the subsequent softening event. During the HSS, 
the compac t core jet is suppressed by over two orders of 
magnitude (|Coriat et al.|[201ll ; iRussell et al.ll201ll ). and re- 
forms in the final stages of the outburst as the source moves 
from a HIMS back to the LHS. 

While many authors have used this model as a bench- 
mark against which to compare their ob s ervat io nal data (e.g. 
Corbel et ail |2004| ; iMiller- Jones et al.l [20051; iFender et all 



20061 ). it was evaluated in detail bv lFender et all ([20091 ), who 



attempted to couple the moment of jet ejection with partic- 
ular X-ray spectral and timing features. While they verified 
the main features of the previous picture, they were unable 
to identify any unique X-ray signature of the jet ejection 
event (the crossing of the jet line), probably due to the sparse 
radio coverage for most of the studied events. These authors 
identified a need for more intensive radio coverage during 
the transitional phases of the outburst, particularly in the 
form of high angular resolution VLBI monitoring. Such re- 
solved radio imaging is crucial in accurately interpreting the 
jet emission, since the peak in the radio light curve may not 
correspond to the exact time when the moving knots were 
eje cted from the core. A c lear example of this was identified 
bv lVermeulen et al.l ([19871 ), who found that the jet emission 
in SS 433 reached a maximum when crossing a 'brightening 
zone' located ~ 250 au downstream from the core. Absorp- 
tion effects, recollimation shocks, or the time delay required 
to form internal shocks could all prevent an accurate as- 
sessment of the disc-jet coupling when relying on radio flux 
density monitoring alone. 

1.1 H1743-322 

First detected durin g an outburst in 1977 
|Kaluzienski fc Holt|[i977l ). H1743- 322 was rediscovered as a 



new source in 2003 by INTEGRAL jRevnivtsev et al.ll2003h . 
only subseque ntly being associated wit h the previously 
known source (Mar kwardt fc Swan kl l2003h . Since that time, 
seven smaller o utbursts of the sour c e have been ob served 
dSwankl [200i; iRupen et all 120051; iKalemci et~aH 



200S 



Kuulkers et alj 120081; iKrimm et all 120091 ; iNakahira et~al] 



20ld ; iKuulkers et all l201ll ), at least one of which failed 
to make a full transit ion through all the canonical states 
([Capitanio et aI.ll2009l V 

From its similarities with the dynamically confirmed 



black hole X-ray binary XTE J1 550-564, the source wa s 
classified as a black hole candidate ([McClintock et al1l2009l ) , 
although no mass function has yet been reported for the sys- 
tem. From its location only 3.2° from the Galactic Centre, 
the source distance w as assumed to be ~ 8kpc. This wa s 
recently confirmed by ISteiner. McClintock fc Reidl (|2012l ). 
who derived a distance of 8.5 ± 0.8 kpc by modelling the 
trajectories of the rad io jets during the 2003 outburst. 
iMcClintock et all ([20091 ) measured quiescent magnitudes for 
the system of 17.1 at K-band and > 24 at i' band, but 
with no firm distance determination and uncertain contri- 
butions from the jet and outer disc, the spectral type of 
the secondary is also unknown. However, from the rela- 
tion between orbita l period and quiescent X-ray luminosity 
([Garcia et al.|[200ll ). and their measure d quiescent luminos - 
ity of ~ 3 x lO^ergs" 1 (0.5-10 keV), Ijonker et all ([20101 ) 
inferred an orbital period of > 10 h. 

From the X- ray dipping behaviou r observed by 
iHoman et all i[2005l ) and I Miller et al.l ([20061 ). the inclination 
angle of the disc in H1743-322 is believed to be relatively 
high, on the order of ~ 80° to our line of sight. ISteiner et all 
(|2012l ) also determined the inclination angle of the radio jets 
to be 75 ± 3°, implying that the disc and jet axes are there- 
fore relatively well aligned. Narrow, variable and blueshifted 
absorption line s from highly-ioniz ed Fe suggest a disk wind 
in this system (|Miller et alj|200~6i ), which was subsequently 
found to be partially state-dep endent, with the wind being 
stronger in the high/soft state (|Blum et al.ll20l51 ). The mass 
outflow rate in the wind was estimated as being on the order 
of 3 x 10 17 gs _1 , assuming a wind velocity of 300 km s -1 . 

The X-ray and radio behaviour of the source during the 
2003 outburst were studied in detail by IMcClintock et all 
([2009h . The event was complicated, with multiple transi- 
tions between the different X-ray spectral states, and sev- 
eral different radio flaring events. Owing to the rapid X-ray 
variability, complicated outburst behaviour and sparse ra- 
dio sampling, no detailed compariso n of the jet-di s c coup ling 
with the phenomenological model of lFender et al.l (|2004h was 
possible. 

Several months after the 2003 outburst, moving radio 
emission was detected at arcsecond scales on either side of 
the central source, s ubsequently discove red to be associated 
with X-ray sources (|Corbel et al.ll2005h . These were inter- 
preted as the interaction of the ejecta from the original out- 
burst with the surrounding environment. While the jets were 
found to have decelerated since the 2003 event, the obser- 
vations placed constraints on the position angle of the jets 
and on their apparent velocity. 

In this paper, we present detailed radio observations 
of the 2009 outburst of H1743-322. The X-ray evolu- 
tio n of the outburst has already been a nalys ed in detail 
by iMotta. Munoz-Darias fc Bellonil (|2010l ) and IChen et al] 
([20101 ), and the source was found to follow the canonical 
path through the HID, with a relatively simple evolution as 
compared to the 2003 event. The dense radio coverage at 
high angular resolution, coupled with the simple, canonical 
behaviour in the X-ray band, allowed us to make a detailed 
study of the disc-jet couplin g in this system, testing the 
model of lFender et al. I (120041 ) and addressing some of the 
issues raised by Fender et al.l (|2009l ) . 
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2 OBSERVATIONS AND DATA REDUCTION 

The outburst was monitored with the Proportional Counter 
Array (PCA) on board the Rossi X-ray Timing Explorer 
(RXTE), and its X -ray spectral and ti ming characteristics 
have been used by iMotta et al.l |2010l ) to distinguish the 
various characteristic accretion states. We monitored the 
outburst in the radio band using the VLA and Very Long 
Baseline Array (VLBA) under the auspices of the Jet Accel- 
eration and Collimation P robe of Transient X-ray binaries 
(JACPOT XRB) program i|Miller- Jones et alj2011af ). track- 
ing the evolution of the radio emission through the entire 
outburst. We also obtained data from the ATCA to improve 
our temporal coverage in the radio band. 

2.1 VLA 

After Swift reported a re-brightening of H 1743-322 in th e 
15-50 keV band from 2009 May 22 (|Krimm et al.l l2009h . 
we made preliminary radio observations with the VLA, 
detecting the source a t a level of 2.2 mjy at 8.4 GHz 
l|Miller- Jones et al.|[2009h . On the basis of this detection, we 
triggered a long-term VLA monitoring campaign (program 
code AM991), tracking the radio emission from the source 
over 24 epochs as it made the transition from the LHS to 
the HSS, and again as it made the reverse transition back 
to the LHS (see Table [1] for more details). The array was in 
its CnB and C configurations during the observations, with 
the configuration change occurring on 2009 June 16. The 
monitoring was carried out at frequencies of 1.4, 4.9 and 
8.4 GHz, with two epochs of observation at 22.4 GHz dur- 
ing the peak of the radio outburst in an attempt to better 
constrain the source spectrum when the source was bright- 
est. We observed in dual-polarization mode with 100 MHz 
of bandwidth per polarization, split equally between two 
intermediate frequency (IF) pairs. We used 3C 286 as a pri- 
mary calibrator, setting the amplitude scale according to 
the coefficients derived at the VLA in 1999 by staff at the 
National Radio Astronomy Observatory, as implemented in 
the 31Dec08 version of the Astronomical Image Processing 
System (MPS; iGreisenl [20031 ). The secondary calibrator was 
J1744-3116, at an angular separation of 1° from H1743-322. 
Observations were carried out in fast-switching mode with a 
4-min cycle time (1 min on the calibrator, 3 min on the target 
source) to minimize down time while maximizing the accu- 
racy of the phase transfer at the low elevation of the sources. 
Data reduction was carried out according to standard pro- 
cedures within AlPS. The observations were typically of du- 
ration 1 h split across all observing frequencies, and were 
simultaneous with monitoring at the VLBA wherever pos- 
sible. On a few occasions, the weather was sufficiently bad 
that the calibrator phase changed by several tens of degrees 
over the course of a 60-s scan. At these epochs, the inter- 
polation of phases to the target source was poor, and since 
H1743-322 was not bright enough for self-calibration we were 
unable to correct for this effect, leading to larger than ex- 
pected errors on the measured source brightnesses. In one 
instance, on 2009 July 4, the weather was sufficiently bad 
(phase changing by > 90° within 60 s) that the target source 
could not be detected. 



2.2 VLBA 

On detection of a radio source with the VLA, we trig- 
gered a campaign of observations with the VLBA under pro- 
gram code BM308, monitoring the source over ten epochs 
throughout its transition from the LHS to the HSS. Once our 
VLA monitoring showed that the radio emission had faded 
below detectable levels for the VLBA, the high-resolution 
monitoring ceased until X-ray monitoring showed the onset 
of the reverse transition back to the LHS, whereupon we 
made a total of six further observations. 

We observed in dual polarization mode at a primary 
observing frequency of 8.4 GHz, although we also observed 
at 15 GHz for the four epochs during the peak of the out- 
burst when the radio emission was brightest. Each epoch 
was of duration 5 h, the maximum time for which the source 
was above the elevation limit of the antennas. For all epochs 
except the last one, our recording rate was 256 Mbps, corre- 
sponding to a total bandwidth of 32 MHz per polarization. 
For the final epoch, when our VLA monitoring showed that 
the source was well into the decaying LHS and very faint, we 
doubled the recording rate and observing bandwidth. In the 
last ten epochs of observation, we spent 30 min at the be- 
ginning and end of the observing run taking data on bright 
extragalactic calibrators distributed over a wide range of el- 
evations across the sky (i.e. geodetic blocks). This enabled 
us to better account for unmodeled clock and tropospheric 
phase errors using the AlPS task DELZN, thereby improving 
the success of the phase transfer. 

Owing to its location only 3.2° away from the Galactic 
Centre, the angular broadening caused by the strong scat- 
tering along the line of sight to the source creates a dearth 
of good, close VLBI calibrators, particularly at lower radio 
frequencies. Based on experience from previous observing 
campaigns on H1743-322 in 2003, 2004 and 2008, we used 
the sources J1748-2907 and J1740-2929 as our primary phase 
reference sources, both believed to be extragalactic on ac- 
count of their brightness temperatures, morp hologies and 
steady fluxes (|Bower. Backer fc Sramekjl2001|). For J1748- 
2907, we used the updated position of iReid fc Brunthalerl 
(|2004 ). Our phase referencing cycle time was 3 minutes, 
spending 2 minutes on H1743-322 and 1 minute on the cal- 
ibrator source during each cycle. Every seventh cycle, we 
substituted a scan on the target source for a scan on the 
check source J1744-3116 from th e third extension t o the 
VLBA Calibrator Survey (VCS-3; [petrov et al.ll2005h , and 
then observed the fringe finder source J1733-1304. Our as- 
sumed calibrator positions are given in Tabled Despite the 
greater calibrator throw to J1748-2907 and J1740-2929, they 
were preferred to J1744-3116 as phase reference sources since 
the latter source is known to be highly scattered at 8.4 GHz, 
reducing the number of antennas for which phase solutions 
could be derived. 

The data were correlated using the VLBA-DiFX 
implementat i on o f the software correlator developed by 
iDeller et all (|2007t ). and subsequently reduced according to 
standard procedures within AlPS. We corrected the cor- 
relator model for the updated Earth orientation parame- 
ters and ionospheric effects, using the system temperatures 
to calibrate the amplitude scale before deriving instrumen- 
tal phase solutions using the fringe finder source. Fringe 
fitting was then performed on the phase reference source 
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i\/r Tri a 


X-ray 


•3 1.4GHz 


•34.9GHz 


-38.4GHz 


^22.4GHz 


"8.4,4.9 


C*8.4,1.4 






(A \ 
W 


state 


(m Jy bm ) 


(nijybm - ) 


(mjy bm - ) 


(m Jy bm - ) 






2009 


May 27 


54978.4 


_ 


_ 


_ 


2.24 ±0.05 


_ 


_ 


_ 


2009 


May 30 


54981.4 


HIMS 


2.03 ±0.30 


2.76 ± 0.10 


2.73 ±0.01 


_ 


-0.02 ±0.09 


0.16 ±0.08 


2009 


Jun 2 


54984.3 


HIMS 


1.80 ±0.30 


2.50 ± 0.06 


2.48 ± 0.04 





-0.01 ± 0.05 


0.18 ±0.09 


2009 


Jun 7 


54989.2 


HSS 





5.7 ± 0.4 


12.8 ±0.3 





1.43 ± 0.13 





2009 


Jun 9 


54991.4 


SIMS to HSS 








2.20 ±0.80 


< 1.41 c 





_ 


2009 


Jun 11 


54993.3 


HSS 


1.71 ±0.30 


0.90 ± 0.08 


0.65 ±0.04 


< l.ll c 


-0.58 ±0.20 


-0.53 ±0.10 


2009 


Jun 14 


54996.2 


HSS 


< 1.95 c 


0.43 ± 0.10 


0.18 ±0.06 





-1.56 ±0.78 





2009 


Jun 15 


54997.3 


HSS 





0.35 ±0.04 


0.38 ±0.06 





-0.18 ±0.35 





2009 


Jun 16 


54998.3 


HSS 





_ 


0.31 ± 0.05 











2009 


Jun 17 


54999.2 


HSS 








0.16 ±0.05 


_ 




_ 


2009 


Jun 18 


55000.3 


HSS 


_ 


_ 


0.17 ± 0.06 


_ 


_ 


_ 


2009 


Jun 19 


55001.3 


HSS 


_ 





0.19 ±0.06 











2009 


Jun 23 


55005.3 


HSS 





_ 


0.34 ± 0.07 


_ 








2009 


Jun 25 


55007.2 


HSS 





_ 


0.35 ±0.06 


_ 







2009 


Jun 27 


55009.3 


HSS 








0.33 ±0.06 










2009 


Jun 29 


55011.2 


HSS 








0.24 ±0.06 











2009 


Jun 30 


55012.2 


HSS 








0.27 ± 0.06 


_ 


_ 


_ 


2009 


Jul 7 


55019.2 


HIMS 


- 


- 


0.59 ±0.06 


- 


- 


- 


2009 


Jul 8 


55020.3 


HIMS 


0.95 ± 0.22 


0.65 ± 0.08 


0.55 ± 0.06 




-0.32 ± 0.30 


-0.31 ± 0.14 


2009 


Jul 9 


55021.3 


HIMS 


0.99 ±0.26 


0.48 ± 0.09 


0.41 ± 0.07 




-0.27 ± 0.47 


-0.49 ±0.18 


2009 


Jul 12 


55024.3 


HIMS 


1.30 ±0.27 


0.48 ± 0.09 


0.34 ± 0.06 




-0.66 ±0.50 


-0.75 ±0.16 


2009 


Jul 14 


55026.2 


HIMS to LHS 


0.71 ± 0.25 


0.62 ± 0.08 


0.59 ±0.07 




-0.09 ±0.32 


-0.11 ±0.20 


2009 


Jul 20 


55032.2 


LHS 


1.25 ±0.27 


0.47 ± 0.08 


0.63 ±0.05 




0.52 ± 0.35 


-0.38 ±0.13 


2009 


Aug 6 


55049.2 


LHS 




< 0.20 c 


0.19 ±0.05 









Table 1. Log of the VLA observations, showing for each epoch the X-ray state of the source, the brightness of H1743-322 at each 
frequency, and the derived radio spectral index a where appropriate. 

a The reported MJD corresponds t o the midpoint of th e observations. 

b X-ray states are those derived by iMotta et al.l ll2010T l . Where the radio data were taken between two X-ray observations in different 
states, we report the relevant transition. 

c In the case of non-detections, we quote the 3<r upper limit to the source brightness. 



Calibrator Right Ascension Declination Calibrator 

throw (°) 

J1748-2907 17 h 48 m 45!6860 -29°07'39 / /404 3.2 
J1740-2929 17 h 40 m 54?5275 -29°29'50'/327 3.0 
J1744-3116 17 h 44 m 23?5784 -31°16'36y2913 1.0 



Table 2. Assumed positions of the phase reference and check 
sources (all in J 2000 coordinates). The last column shows the 
angular separation between the calibrator and H1743-322. 

(the brighter source J 1748-2907 where available, otherwise 
J1740-2929), which was subjected to iterative imaging and 
self-calibration. The final image of the phase reference cali- 
brator was then used as a model for bandpass calibration 
before transferring bandpass, amplitude and phase solu- 
tions to H1743-322. However, the large calibrator throw (Ta- 
ble [2]) coupled with the low elevation of the sources meant 
that unsolved phase gradients compromised the astromet- 
ric accuracy of the observations, with the positions of both 
H1743-322 and the check source being found to vary between 
epochs. Since this had a detrimental effect on our ability to 
measure the proper motions of the jet ejecta, we tested a hy- 
brid approach which was found to give a much more stable 
position of the core radio emission from H1743-322 during 
the LHS in the initial phase of the outburst. We first solved 



for the short timescale phase variations using the original 
phase reference source. Having transferred the solutions to 
H1743-322, we then re-derived the phase solutions using the 
check source J1744-3116 as the phase calibrator. Owing to 
the infrequent observations of this source, this imposed a 
longer-timescale phase variation on the solutions. The im- 
portant effect was to re-reference the astrometry to provide a 
more stable target position relative to the assumed position 
of the closer calibrator J1744-3116. 

Owing to the strong scattering along the line of sight, we 
were restricted to baselines shorter than 25 MA at 8.4 GHz, 
implying that we could only derive meaningful information 
from baselines between the five south-western antennas (FD, 
KP, LA, PT, OV). The loss of half our antennas both re- 
duced our sensitivity by a factor of 2 relative to that pre- 
dicted for the full array, and degraded our angular resolution 
and Mf-coverage. Nevertheless, we were able to image the 
data, using natural weighting to maximize our sensitivity. 
The results are summarized in Table 

2.3 ATCA 

Between 2009 May 27 (MJD 54978.8) and 2009 August 7 
(MJD 55050.34), we performed a total of 8 observations 
of H1743-322 with the Australia Telescope Compact Array 
(ATCA). The observations were carried out using the new 
Compact Array Broadband Backend (CABB; I Wilson et al.l 
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Date 




i\/r Tri a 


X-ray 


Frequency 


Bandwidth 


Phase reference 


Geodetic 


ril(4o-ozz nUX 






1 . \\ 

w 


state 


[yjiiz ) 


^IVlriz ) 


calibrator 


DIOCK I 


density c (mjy) 


2009 


May 28 


54979.35 ± 0.10 


_ 


8.4 


32 


J1748-2907, J1740-2929 


N 


2.0 ±0.2 


2009 


May 30 


54981.35 ± 0.10 


HIMS 


8.4 


32 


J1748-2907, J1740-2929 


N 


2.1 ± 0.3 


2009 


Jun 2 


54984.35 ± 0.10 


HIMS 


8.1 


32 


J1748-2907, J1740-2929 


N 


1.8 ±0.3 


2009 


Jun 5 


54987.33 ± 0.10 


SIMS 


8.4 


32 


J1748-2907, J1740-2929 


N 


5.7 ± 0.3 






54987.33 ± 0.10 


SIMS 


15.0 


32 


J1748-2907, J1740-2929 


N 


2.2 ±0.4 


2009 


Jun 6 d 


54988.33 ± 0.10 


HSS 


8.4 


32 


J1740-2929 


N 


6.4 ± 0.3 








HSS 










7.1 ± 0.3 






54988.34 ± 0.09 


HSS 


15.0 


32 


J1748-2907 


N 


2.9 ±0.5 


2009 


Jun 8 


54990.33 ± 0.10 


SIMS 


8.4 


32 


J1740-2929 


N 


< 1.0 






54990.34 ± 0.09 


SIMS 


15.0 


32 


J1740-2929 


N 


< 1.7 


2009 


Jun 9 


54991.33 ± 0.08 


HSS 


8.4 


32 


J1740-2929 


Y 


< 1.1 






54991.33 ± 0.09 


SIMS to HSS 


15.0 


32 


J1740-2929 


Y 


< 1.1 


2009 


Jun 11 


54993.33 ± 0.08 


HSS 


8.4 


32 


J1740-2929 


Y 


< 0.8 


2009 


Jun 13 


54995.33 ± 0.08 


HSS 


8.4 


32 


J1740-2929 


Y 


< 1.0 


2009 


Jun 15 


54997.31 ± 0.08 


HSS 


8.4 


32 


J1740-2929 


Y 


< 0.8 


2009 


Jul 9 


55021.23 ± 0.09 


HIMS 


8. 1 


32 


J1748-2907 


Y 


< 0.8 


2009 


Jul 11 


55023.25 ± 0.08 


HIMS 


8.4 


32 


J1748-2907 


Y 


< 0.8 


2009 


Jul 12 


55024.25 ± 0.08 


HIMS 


8.4 


32 


J1748-2907 


Y 


< 0.8 


2009 


Jul 18 


55030.23 ± 0.08 


LHS 


8.4 


32 


J1748-2907 


Y 


< 0.8 


2009 


Jul 20 


55032.06 ± 0.06 


LHS 


8.4 


32 


J1748-2907 


Y 


< 0.9 


2009 


Aug 6 


55049.18 ± 0.09 


LHS 


8.4 


64 


J1748-2907 


Y 


< 0.4 



Table 3. Log of the VLBA observations, showing for each epoch the X-ray state of the source, the observing frequencies and bandwidths, 
phase calibrators, and whether a geodetic block was used. 

a The reported MJD corresponds to the midpoint of the observations. Error bars denote the duration of the observations cither 
side of the midpoint, such that the total duration corr esponds to twice the reported time uncertainty. 
b X-ray states are those derived by iMotta et al ] l|20ld l. 

c Where the source was not detected, 5<r upper limits are given in all cases. This shows the confidence with which wc can rule out the 
detection of jet ejecta at an unknown position. However, detection of the core can be ruled out at the 3cr level, i.e. 40 per cent lower 
than the tabulated figures, since the core position is known. 

d At this epoch, 2 components were detected at 8.4 GHz, on either side of the core position (see Fig. [TJ. Flux densities for both 
components are reported. At 15.0 GHz, only one component was detected, whose position was consistent with the western component 
in the 8.4 GHz image. 



l201ll ). This upgrade has provided a new broadband backend 
system for the ATCA, increasing the maximum bandwidth 
from 128 to 2048 MHz. Each observation was conducted si- 
multaneously in two different frequency bands, with typical 
central frequencies of 5.5 GHz and 9 GHz. Various array 
configurations were used during our observing campaign. On 
May 27 and July 17 the array was in the compact, hybrid 
(north-south and east-west baselines) configurations H214 
and H75 respectively. The more extended configurations 6A 
and 6D were used for the remaining 6 observations. 

The ATCA has orthogonal linearly polarized feeds and 
full Stokes parameters (I, Q, U, V) were recorded at each 
frequency. We used PKS 1934-638 for absolute flux and 
bandpass calibration, and PMN 1729-37 to calibrate the 
antenna gains and phases as a function of time. We deter- 
mined the polarization leakages using either the primary 
or the secondary calibrator, depending on the parallactic 
angle coverage of the secondary. Imagin g was carried out 
using a multi-frequency clean algorithm |Sault fc Wieringal 
1 1994! ) to take into account spectral variations of the source 
flux density across the 2 GHz bandwidth. The editing, cal- 
ibration, Fourier transformation, deconvolution and image 
analysis were carried out with the Multichannel Image Re- 
construction, Image Analysis and Display (miriad) software 



l|Sault. Teuben fc Wrightl I1995T I. The results of the ATCA 
observations are given in Table 0] 

2.4 RXTE 

We analysed all publicly available observations of the 2003 
and 2009 outbursts of H1743-322 in the RXTE archive. 
The data were reduced using the HEASOFT software pack- 
age v6.8, following the standard steps described in the 
RXTE cookbookfl We constructed Hardness Intensity Dia- 
grams (HIDs) from the Proportional Counter Array (PCA; 
IJahoda et al.l 120061 ) data. We produced background sub- 
tracted light curves binned at 16 s using PCA Standard 2 
mode data from Proportional Counter Unit (PCU) 2 (all lay- 
ers). The light curves were divided into three energy bands, 
2.1-4.9, 8.6-18.3 and 2.1-18.3 keV (according to the corre- 
spondence between fixed energy channels and energy valid 
for PCA gain epoch 5). We defined the hardness ratio (HR) 
as the ratio of count rate in the bands 8.6-18.3 keV and 
2.1-4.9 keV and the intensity as the count rate in the en- 
ergy band 2.1-18.3 keV. 

All the observations we used fall in PCA gain epoch 

1 http://hcasarc.gsfc.nasa.gov/docs/xte/data_analysis.html 
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Figure 1. 8.4-GHz VLBA images of H1743-322. Contour levels are at ±(v / 2) n times the rms noise level indicated in the top left corner of 
each image, where n = —3, 3, 4, 5... The grey scale shows the image flux density in units of mjy beam -1 , except for the stacked image of 
2009 July 9-18, for which the units are /xjybeam -1 . The white cross marks the weighted mean core position derived in Scction l3.3.1l We 
see an unresolved core jet in the first three epochs, followed by a two-sided ejection event. The positional offset of the possible northwestern 
source in the image of 2009 July 9-18, together with the marginal nature of the detection, implies that we cannot convincingly associate 
this source with the compact core jet. 



5 (from 2000 May 13 until the present). However, (small) 
continuous gain changes are reported within one epoch. To 
correct for this effect, which could affect our comparison 
of the 2003 and 2009 outbursts, we estimated the varia- 
tion of the HR of the Crab between 2003 and 2009. Be- 
cause the spectrum of the Crab is supposed to be constant, 
any changes observed in the hardness of the Crab are most 
probably caused by changes in the instrumental response. 
We selected RXTE observations covering the same period 
as our observations of H1743-322. We estimated an increase 
of ~ 5.7 per cent in the average HR of the Crab from 2003 
to 2009. Consequently, we divided all the HRs calculated for 
H1743-322 in 2009 by a factor 1.057 to match the 2003 data. 
The changes in the PCA gain, although small, might also af- 
fect the source count rate. On average, the Crab count rate 
in the 2.1-18.3 keV decreased by ~ 4.6 per cent from 2003 
to 2009. We thus corrected the intensities of HI 743-322 to 
account for this change. 

To determine the spectral state of the source at the time 
of each RXTE ob servation, we used the classifications of 
iMotta et ail (I^Old ). These were based on timing properties 



(the fractional rms variability, the presence of QPOs, and the 
overall shape of the power density spectra), hardness ratios, 
and also on the relative contributions of disc and power- 
law components to the X-ray spectra (see for example their 
fig. 5). This latter diagnostic was important in distinguishing 
the SIMS (with significant power-law contribution) from the 
HSS (dominated by the thermal disc emission), since the 
fractional rms variability and hardness ratios could be very 
similar in these two states (see Fig. [2}. 

3 RESULTS 

3.1 Radio light curves and spectra 

The radio light curves and the evolution with time of the 
radio spectral index a (defined as S v oc v a ) are shown in 
Fig. [2] Also plotted for comparison purposes are the X-ray 
count rate, hardness and integrated fractional rms variabil- 
ity from RXTE/PCA monitoring, togeth er with the X-ray 
state classification s oflMotta et all (|2010T ). As expected (e.g 
I Fender et al. | |2004 120091 )7 changes in the radio emission are 
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Date 


MJD a 


A-ray 


i_ 1LUUL11L \ 


J_ 1 U..A. LlCllGlL y 








statp^ 


(GHz) 


( tyi l\z m 1 ^ C 

l 111U y L/OCLJLIX } 




2009 May 27 


54978.8 


— 


5.50 


2.40 ±0.07 


-0.28 ± 0.08 






— 


8.62 


2.12 ± 0.05 


— 


2009 Jun 4 


54986.6 


HIMS to SIMS 


5.50 


0.26 ±0.02 


-0.73 ± 0.31 






— 


8.80 


0.18 ±0.03 


— 


2009 Jun 5 


54987.6 


SIMS to HSS 


5.50 


23.82 ± 0.14 


0.03 ±0.02 






— 


9.00 


24.20 ±0.11 


— 


2009 Jun 6 


54988.8 


HSS 


5.50 


24.37 ± 0.09 


0.15 ±0.01 






— 


9.00 


26.19 ± 0.16 


— 


2009 Jun 17 


54999.8 


HSS 


5.50 


0.23 ±0.03 


-0.66 ± 0.45 






- 


9.00 


0.17 ± 0.03 


- 


2009 Jun 20 


55002.6 


HSS 


5.50 


0.20 ± 0.02 


—0.67 ± 0.48 








9.00 


0.14 ±0.03 




2009 Jul 13 


55025.4 


HIMS to LHS 


5.50 


< 0.16 










9.00 


< 0.30 




2009 Aug 7 


55050.8 


LHS 


5.50 


0.10 ±0.02 










9.00 


< 0.11 





Table 4. Log of the ATCA observations, showing for each epoch the X-ray state, the measured flux densities, and the derived radio 
spectral index of H1743-322. 



a The reported MJD corresponds t o the midpoint of th e observations. 

b X-ray states are those derived by iMotta et al ] l|20ld l. Where the radio data were taken between two X-ray observations in different 
states, we report the relevant transition. 

c In the case of non-detections, we quote the 3<r upper limit to the source brightness. 



more closely coupled with X-ray spectral and timing infor- 
mation than with the overall X-ray light curve. 

By the time we triggered our observing campaign, the 
source was already in the HIMS, with a steady radio flux 
density of ~ 2 mjy and a flat radio spectrum. This per- 
sisted until the end of the HIMS, after which the ATCA de- 
tected a quenching of the radio emission and a steepening of 
the radio spectrum to a = -0.76 ± 0.32 on MJD 54986.6. 
Quenching of the radio emission prior to a major ejec- 
tion event is a well- known phenomenon in Cygnus X-3 
(|Waltman et al.lll994T), and has also been observed in XTE 
J1859±226 (|Brocksopp et al.ll2002T ). associated with a brief 
steepening of the radio spectrum. It is likely that this ob- 
servation corresponds to the phase of instability i n the jet 
prior to the main outburst noted by I Fender et al.l (|2004l ). 

Within 0.8 d of the ATCA detection of quenched ra- 
dio emission, the VLBA detected the onset of the radio 
flare, on MJD 54987.3. This observation was almost simul- 
taneous with the first X-ray detection of the SIMS. The 
peak of the radio outburst, as detected with the ATCA, 
reached 26.2 mjy on MJD 54988.8, and the radio brightness 
as measured by all instruments then decreased for the next 
14 d as the source moved into the HSS. At the peak of the 
flare, the radio spectrum was slightly inverted, consistent 
with expanding synchrotron plasma that was self-absorbed 
at lower frequencies, with a turnover in the observed ATCA 
frequency range. 

While RXTE detected an excursion from the HSS back 
into the SIMS on MJD 54990.3, the radio sampling was too 
coarse to detect a temporary rebrightening at this time. The 
VLBA did not detect the source on MJD 54990.3, and the 
measured VLA radio brightness on MJD 54991.4 was con- 
sistent with the decay from the original flare. The only ev- 
idence for a second ejection event is the inverted spectrum 



between 4.9 and 8.5 GHz on MJD 54989.2 (a = 1.4 ± 0.2). 
Given that these data were taken after the peak of the main 
flare as measured by ATCA, this suggests the early stages 
of a second ejection event. 

Following this VLA detection of an inverted spectrum, 
the radio spectra were steep during the HSS. The source 
was detected throughout this phase of the outburst, with a 
slight increase in the brightness at 8.5 GHz on MJD 55005.3, 
following which the radio emission was relatively stable at a 
level of ~ 0.3 mjy. The unified model of lFender et~aH (|2004T l 
suggests that the compact core jet should be suppressed in 
this phase, so one likely explanation for the low-level emis- 
sion observed at late times in the HSS is an ongoing interac- 
tion between the jet and the interstellar me dium through 
which it was propagating (as suggested by ICorbel et al.l 
|2004| . to explain the residual radio emission detected in the 
HSS of XTE J1650-500). However, the source was not de- 
tected with the VLBA during this period, and the relatively 
compact configurations of the VLA and the ATCA did not 
provide sufficient resolution to distinguish between any such 
external shocks and a compact core jet. Thus in the absence 
of either spectral or imaging information, we cannot verify 
this hypothesis. 

After the reverse transition to the HIMS, the radio 
brightness showed a second slight increase, to ~ 0.6 mjy, 
with a spectrum consistent with being flat. We interpret 
this as the reignition of the compact core jet. As the X-ray 
emission decayed late in the outburst, the radio brightness 
also decreased, reaching 0.1 mjy by MJD 55050.3. 



3.2 Radio imaging 

H1743-322 was detected by the VLBA on 5 occasions during 
the outburst (Table[3]). Three of these were during the HIMS 
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h. Hi-. 




1 




Time (MJD-54975) 



3.2.1 No evidence for diffuse, extended emission 

In the initial LHS, the flux density measured with the VLB A 
was at a level consistent (at the 2cr-level) with that ob- 
served during the quasi-simultaneous VLA and ATCA ob- 
servations. However, the two VLBA observations during the 
rise and at the peak of the flare recovered only 25-50 per 
cent of the integrated flux density measured a few hours 
later by the ATCA. Although this could indicate excess dif- 
fuse emission that was resolved out by the VLBA, it is likely 
that the majority of the discrepancy can be attributed to ei- 
ther timing effects during the rise phase of the flare, or to 
atmospheric phase decorrelation. The VLBA observation on 
MJD 54987.33 ± 0.10 was taken during the rise phase of the 
flare, 6 hours prior to the closest ATCA observation. The 
ATCA data (Table Q} show that the 9-GHz radio emission 
increased from 0.2 to 24.2 mjy over the 24-hour period from 
MJD 54986.6 to 54987.6, so a 6-hour time offset could cor- 
respond to a significant change in flux density. 

To test the effects of phase decorrelation, we com- 
pared the flux density recovered from the calibrator sources 
J1748-2907 and J1740-2929 after applying the phase cor- 
rections derived using J1744-3116, with that derived from 
self-calibration of those same two sources. This showed that 
phase decorrelation could reduce the recovered flux density 
by up to 50 per cent on MJD 54987.3, and by up to 25 per 
cent on MJD 54988.3 (when the atmosphere was presumably 
somewhat more stable). Taking into account the effects of 
both time offsets and atmospheric phase decorrelation, we 
conclude that although we cannot rule out the existence of 
diffuse emission that is resolved out by the VLBA, such a 
component is not required to reconcile the measured VLBA 
fluxes with those seen at the ATCA. 



Figure 2. Evolution of the 2009 outburst of H1743-322. (a) 2.1- 
18.3 keV count rate as measured by RXTE/PCA. (b) X-ray hard- 
ness, defined as the ratio of the count rate in the 8.6-18.3 keV 
band to that in the 2.1-4.9 keV band, (c) X-ray timing informa- 
tion. The plot shows the integrated fractional rms variability in 
the 0.1—64 Hz band, with symbols rep resenting notable fe atures in 
the power spectrum. Data taken from lMotta et all l|2010h . (d) ra- 
dio light curves, (e) radio spectral indices. The grey shaded bands 
show the X-ray sta tes determined for each RXTE observation by 
iMotta et al ] l|2010l l. as labelled in panel (d). The yellow band 
(bordered by the vertical dashed lines) shows our inferred ejec- 
tion date for the relativistically- moving knots (see Section l3.3.2|l . 



prior to the radio flare, and the following two were during 
the rise phase and at the peak of the flare. The resulting im- 
ages, together with the stacked image from all data during 
the outburst decay, are shown in Fig. [T] In the first three 
observations during the HIMS, we detected only the com- 
pact core jet, which was unresolved down to a beam size 
of 14.8 x 5.6 mas in PA. 6.6° E of N. In contrast, at the 
peak of the flare, we observed two distinct components, sug- 
gesting a two-sided ejection event. We note that the weak 
component seen at the eastern edge of the VLBA image from 
2009 June 5 (with relative co-ordinates (65,-23) mas in R.A. 
and Dec, respectively) has a significance of < 5<r. Given 
the marginal nature of this detection, and since its orienta- 
tion is not aligned with the observed axis of the bright jet 
components, we do not believe it to be real. 



3.2.2 No evidence for free-free absorption in the inner 
regions 

Given the hi gh inclination of th e system and the presence of 
a disk wind (|Miller et alj|2006l ). we should address whether 
there is any evidence for free-free absorption in the inner 
radio-e mitting regions, a s has been found in the case of 
SS433 (|Paragi et al.ll 19991 ). Since the outflow rate estimated 
for the wind in H1743-322 (|Miller et al.ll2006l ) is several or- 
ders of magnitude lower than that in the equatorial wind of 
SS433 (|Blundell et al.ll200ll ), there should be significantly 
less free-free absorbing material along our line of sight to 
the inner jets. This is supported by the radio spectrum, 
which is only ever highly inverted (as would be expected 
for free-free absorption) during the rise phase of the flare 
(Fig. [2]), when an inverted spectrum can be attributed to 
low-frequency synchrotron self-absorption in the expanding 
ejecta. Finally, we always detect the core of the system with 
our VLBA imaging during the HIMS, at a position interme- 
diate between the eastern and western jet ejecta seen later in 
the outburst (Fig. [TJ . Thus there is no reason to believe that 
there is significant free-free absorption in the inner regions 
of the system. 
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Figure 3. Fitted positions of all detected VLBA components 
from H1743-322, relative to the weighted mean core position of 
17 h 46 m 15!59637, -32° 14' 00'.' 860 (J 2000). East is to the left and 
north is up. Each point is labelled with the MJD of the observa- 
tion. Note that the positions of both components found for the 
final observation (MJD 54988) arc shown. 



MJD-54975 

Figure 4. Angular separation from core position in R.A. and 
Dec. as a function of time. We also show a fit to the approaching 
component for epochs MJD 54987 and 54988, using this to derive 
an ejection date and hence fit the proper motion of the receding 
component. 



3.3 Astrometry 

3.3.1 Core position 

The early triggering of the observations and consequent de- 
tection of the compact core jet provided a point of reference 
for the subsequent observations during the radio outburst. 
We measured the source position, relative to the assumed 
position of the calibrator source J1744-3116 as given in Ta- 
ble [2] The three observations from the HIMS showed a very 
consistent position in R.A. but a scatter of ~ 5 mas in Dec. 
(Fig. O. The calibrator throws to both the short-timescale 
phase reference calibrators J1748-2907 and J1740-2929 and 
to the final calibrator source J1744-3116 are all significantly 
greater in Dec. than in R.A.. Therefore the phase gradients 
with elevation expected for such a low-declination source will 
have a greater effect on the Dec. component of position than 
on the R.A. component, so such a discrepancy in the scatter 
is not surprising. A weighted mean of the three measured 
core positions from the HIMS provides our best estimate of 
the true position of H1743-322, 

RA = 17 h 46 m 15! 59637(6) 
Dec. = -32° 14'00" 860(2) (J2000), 

where the assumed errors are purely statistical. To this we 
must add a systematic error of 2.7 mas in R.A. and 8.0 mas in 
Dec. from combining the quoted uncertainty on the calibra- 
tor positiorQ and the uncertainty from the calibrator throw 
to J1744-3116 (using a linear extra polation of the results 
of iPradel. Chariot fc Lestradd 12006] to the low declination 
of the target source). We note that our quoted position is 
measured relative to the final calibrator source, J1744-3116, 
whose assumed position is given in Table [2] 

3.3.2 Ejection event 

While there is some scatter in the Dec. component of the 
source position during the HIMS, the R.A. is constant to 

2 http:/ /www. vlba.nrao.edu/astro/calib/ 



within 1.5 mas. After the transition to the SIMS and the 
subsequent HSS, the measured source positions at both 8.4 
and 15 GHz on MJD 54987 and 54988 show evidence of com- 
ponents moving away from the X-ray binary (Fig. On 
MJD 54988, we see two radio-emitting components, aligned 
along an axis oriented 87.2 ± 3.6° E of N. This is fully con- 
sistent with the jet axis (89.0 ± 1.5° for th e eastern jet, 
— 91. 7 ± 1.8° for the western jet) determined bv lCorbel et alj 
(|2005l ). suggesting that the jet axis has not precessed since 
2003. Further, their measurement of larger angular separa- 
tions for the eastern component is in agreement with our 
VLBA observations (the eastern component appearing first, 
on 2009 June 5), despite the fact that their jet knots had 
decelerated by the time the radio emission was observed. 
This suggests that the interstellar medium surrounding the 
source is not sufficiently inhomogeneous that the approach- 
ing (eastern) ejecta were decelerated more th an the reced- 
ing (w estern) component in the observations of lCorbel et al.l 
(|2005h . 

The angular separation of the two components from 
the weighted mean core position is shown in Fig. U As- 
suming that we are identifying a single ballistically-moving 
component on each side of the core (likely, given that only 
one day separates the two observations), a linear fit to the 
displacement along the jet axis of the approaching (east- 
ern) component on MJD 54987 and 54988 gives a proper 
motion /i apP = 3.7 ± 0.7masd _1 (an apparent velocity of 
"app = 0.18c at 8.5 kpc). This fit determines the ejection 
date as i e j =MJD 54984.9±o g. While this ejection date was 
derived from only two measurements of the ejecta positions, 
the agreement found between the positions measured in the 
15 GHz and 8.4 GHz VLBA images gives us confidence in 
our estimate. With only two data points, there may be ad- 
ditional systematic uncertainties affecting the determination 
of the ejection date, although their magnitude is not possi- 
ble to quantify from our data. While we cannot definitively 
rule out d eceleration of the ejecta on such small angular 
scales (e.g. lYang et alj|2010l ). we consider it unlikely, since 
the ejecta from the 2003 outburst decelerated much further 
downstream. The continued jet activity of H1743-322 since 
that time (Section ll.lj l makes it unlikely that the density 
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of the interstellar medium close to the central binary could 
have been substantially increased in the interim. 

Assuming therefore that our derived ejection date is 
correct, then a ballistically-moving receding (western) com- 
ponent would have a proper motion /i rcc = 3.3 ±0.8 mas d~ 



0.16c at 8.5 kpc). The proper motions of both ap- 



proaching and receding jet k nots are therefore si gnificantly 
lower than those derived by ICorbel et al.l (|2005t) on larger 
angular scales (21.2 ± 1.4 mas d -1 for the approaching jet, 
13.3±0.6 mas d _1 for the receding jet, although the jets ha d 
likely decelerated by this point; see also lHao fe Zhang|2009l ). 
Either the jet speed varies between outbursts (possibly cor- 
related with jet power; the 2003 outburst was significantly 
brighter in the radio ba nd than the one presented here; see 
iMcClintock et al.l l2009h . the ballistic motion model is not 
applicable to our data, or we are seeing different compo- 
nents in the eastern jet in the two images. 



4 X-RAY DIAGNOSTICS OF THE EJECTION 
EVENT 

A major goal in the study of the coupling between inflow and 
outflow in X-ray binaries is to relate the observed changes in 
the radio emission (probing the relativistic jets) to those seen 
in the X-ray band (probing the inner regions of the accretion 
flow, and possibl y the base of the j ets). The so-called 'jet 
line' proposed by Fender et alj (|2004T l to mark the transition 
from steady, partially self-absorbed, compact core jets to 
bright, transient, relativistically moving ejecta, is a useful 
way to visua l ise th is transition. However, as discussed by 
iFender et al](|2009r ). the reality is likely to be more complex. 
Nevertheless, as we now discuss, even if the exact sequence 
of events has not been uniquely identified to date, certain 
characteristic X-ray spectral and timing indicators appear 
to be closely coupled with the radio ejection event. 



4.1 Spectral indicators 

The evolution of black hole X-ray binary outbursts is often 
visualized in an X-ray hardness- intensity diagram (HID). 
This shows how X-ray intensity variations couple with X- 
ray spectral shape, with harder spectra being dominated 
by the power-law component and softer spectra by black- 
body disc emission. As shown in Fig. H1743-322 under- 
went a fairly typical evolution thro ugh the HID (com pare 
with the schematic model in fig. 7 of lFender et al.|[2004l and 
th e observed behavio ur of the 14 outbursts plotted in fig. 1 
of lFender et~afll2009h • We have overplotted the radio detec- 
tions, showing the variation of the radio brightness as the 
source evolved through the outburst. There is clear evidence 
for radio flaring as the source moves left across the top of 
the diagram from the HIMS to the HSS. Also noteworthy 
are the position of the radio quench phase and the radio 
rebrightening as the source moves back to a harder state 
towards the end of the outburst. 



4.2 Variability indicators 

The transition from the HIMS to the SIMS is defined 
by X-ray timing information, namely a sharp decrease in 
the integrated fractional rms variability and the possible 



300 



100 




-0.8 



0.3 0.5 
Hardness (8.6-18.3 keV)/(2. 1-4.9 keV) 



-0.4 0.0 0.4 0.8 

log(Radio flux density (mjy)) 



1.0 



Figure 5. Hardness-intensity diagram for the 2009 outburst of 
H1743-322. Circles indicate radio detections, with the size and 
shade of the circles indicating the measured radio flux density. 
The location in the HID at the time of each radio observation 
has been interpolated from RXTE PCA observations within 2 d. 
The large star represents the moment of ejection as estimated in 
Section l3.3.2l with the smaller stars showing the extent of the lc 
error bars on this date. Arrows show how the source evolves with 
time. The source follows a fairly canonical path through the HID. 



appearance of Type B QPOs (jBelloni et al.l 120051 ). Using 
case studies of the three sources XTE J155 0-564, XTE 
J1859+226 and GX339-4, IFender et ail (|200gh considered 
whether the sharp changes in X-ray power spectra seen dur- 
ing the transition from the HIMS to the SIMS could be 
used to determine the exact mo ment of the ma j or eje ction 
event, as previously propo s ed bylFerreira et all (l2006h and 
iKlein-Wolt fc van der Krisl (120081 ). However, IFender et all 
(|2009l ) deemed the available radio coverage insufficient to 
associate the exact moment of ejection with any particular 
X-ray spectral or timing characteristics. Sparse time sam- 
pling, coupled with the (unknown) intrinsic delay between 
ejection and radio flaring caused by both optical depth ef- 
fects and the time taken for the development of the hypoth- 
esised internal shocks, made it impossible to determine the 
moment of ejection and hence search for any corresponding 
X-ray signature. With our dense radio coverage, the core lo- 
cation provided by the VLBI detection in the early phase of 
the outburst (Section 13.3. 1 p . and the fitted proper motion 
of the ejecta (Section I3.3.2[) , we finally have a data set in 
which it is possible to probe this connection directly. 

Fig. [6] shows how the X-ray variability properties vary 
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Figure 6. Hardness-rms variability diagram (HRD) for the 2009 
outburst of H1743-322. (a) Measured X-ray data with error bars, 
with symbols re presenting th e featu res in the power spectrum. 
Data taken from lMotta et aL While we note their use of 

different bands for calculating the hardness ratio from those de- 
fined in Section l2.4l this should not affect the overall picture, (b) 
Variation of radio flux density with position in HRD. X-ray spec- 
tral and timing information at the time of each radio observation 
have been interpolated from RXTE PCA observations within 2 d. 
The size and shade of the circles indicate the measured radio flux 
density. The large star represents the moment of ejection as esti- 
mated in Section 13.3.21 with the smaller stars showing the extent 
of the lcr error bars on this date. Arrows show how the source 
moves along the grey trace through the diagram with time. The 
moment of ejection appears to coincide well with the drop in rms 
variability. 



with spectral hardness during the 2009 outburst of HI 743- 
322, and also how the radio emission relates to these diagnos- 
tic X-ray characteristics. Our inferred ejection event occurs 
as the fractional rms variability begins to drop off sharply, 
coinciding with the disappearance of Type C QPOs. Unfor- 
tunately, this phase of the outburst coincides with a 3 d gap 
in the RXTE coverage. Although we do not analyse Swift 
data in this paper, we note that there were no Swift ob- 
servations during the gap in RXTE coverage, so we cannot 
constrain the simultaneity of the ejection event and this X- 
ray signature more accurately. The radio quenching occurs 
as the X-ray variability decreases, but the time taken for 
the development of shocks and their evolution to become 
optically thin in the GHz regime means that the peak radio 
emission is seen only in the zone of reduced variability. 



5 JET PROPERTIES 
5.1 Jet speeds 

Special relativistic effects imply that the ratio of angular 
separations of approaching and receding components from 
the core at any given epoch is 



1 + 13 cos 9 

i - p cos e ' 



(1) 



where /3c is the jet s peed and 9 is the inclinatio n angle to 
the line of sight (e.g. iMirabel fc Rodrfguez1 [l999). From the 
epoch in which we detect both approaching and receding 



components, our best core position implies Z ap p/^rec = 1.12± 
0.13, which gives /3cos(9 = 0.058 ±0.05 6. Combined with th e 
inclination angle of 75 ± 3° derived bv lSteiner et al.1 (|2012f ). 
this implies 0.19 < P < 0.28. 

As noted in Section ^. 3. 21 the proper motions measured 
during the 2009 outburst are significantly lower than those 
measured (albeit on larger scales) in the 2003 outburs t, lead- 
ing to lower derived jet speeds. ICorbel et al.l (|2005l ) found 
a jet speed P = 0.79 (with 6 = 73°) for the 2003 outburst, 
assuming a distance of 8 kpc and ballistic motion of the ap- 
proaching and receding ejecta. Even in the case of deceler- 
ating jet knots, they were able to constrain P > 0.57 ± 0.05. 
Since acceleration of the ejecta during propagation is un- 
likely, and since the consistent position angle of the jet axis 
on the plane of the sky argues against a precessing jet, it 
seems that the jet velocity was very different in these two 
outbursts. Since few X-ray binaries have been observed with 
VLBI during multiple ejection events, we do not have good 
statistics on whether jet speeds are likely to vary between 
different events in a single source, although plausible evi- 
dence for a variable jet speed in GRS 1915+105 has been 
found (|Miller- Jones et alj|2007l ). 



5.2 Suppression of the core jets in the soft state 

The compact core jets are known to b e suppressed in th e 
HSS of H1743-322, by a factor of ~ 700 (|Coriat et alj|201ll ). 
However, our VLA observations during the HSS show low- 
level emission a t 8.4 GHz, w i th a flux density of up to 
0.35 ± 0.06 mjy. ICorbel et ail (|2004h put forward two pos- 
sible explanations to explain similar residual emission de- 
tected in the HSS of XTE J1650-500; namely a residual com- 
pact core jet, and interactions of the ejecta from the outburst 
with the surrounding interstellar medi um. A compact core 
jet should have a flat radio spectrum l|Blandford fc Konigll 
ll979l ; lFenderll200ll) . whereas we would expect optically-thin 
synchrotron emission from a downstre am external shock , 
givin g rise to a steep spectrum (e.g. ICorbel et al. I 120021 , 
120051 ). Unfortunately, we do not have any spectral infor- 
mation in this phase to distinguish between these two pos- 
sibilities, since the source was sufficiently faint that we 
could not split the observing time between multiple frequen- 
cies. However, we can stack our VLBA observations from 
this phase (MJD 54991-54997) to place a 3ct upper limit of 
0.27 mjy beam -1 at our derived core position. Given the un- 
certainties on the HSS flux densities measured by the VLA, 
this is insufficient to conclusively rule out a low-level com- 
pact jet, but, together wit h the known radio suppression 
in the HSS of this system (|Coriat et al.|[201ll ). appears to 
favour the downstream external shock scenario. 



6 THE REVERSE TRANSITION 
6.1 Jet reactivation 

Owing in part to the speed of the X-ray evolution, there 
is relatively little existing radio data constraining the re- 
activation of the jet during the reverse transition from the 
HSS through the HIMS to the LHS. From a c omparison of 
data from two different outbursts of GX339-4, iFender et al.l 
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|2009l ) found radio upper limits during the reverse transi- 
tion at the same hardness where the peak of the flare was 
observed during the original hard to soft transition, imply- 
ing that either the jet line was not vertical or that it was 
not consistent between different outbursts (see Section [7|). 

The reactivation of the jet has been reported from in- 
frare d /optical monitoring of three black hole X-ray bina- 
ries Jjain et al.ll200ll; iBuxton fc Bailvnl [20o3 ; ICoriat et all 
120091 ; iRussell et alj bOlof T Here, the higher energy emis- 
sion from the compact core jet originates in a region close 
to the black hole, within a light second (e.g. ICasella et al.l 
l20ld : iGandhi et aT]l2010l ). and cannot be caused by collid- 
ing ejecta launched previously (assuming a single power law 
with q = —0.6, an extrapolation of the brightest optically 
thin radio ejecta of Hf 743-322 would be < O.f mjy at near- 
infrared (N1R) frequencies; much fainter than the observed 
NIR jet emission in other black hole X-ray binaries at similar 
distances). In all cases, the optical/NIR jet became visible 
only after the source had fully completed the transition to 
the LHS. The radio emission, at least in one of these systems, 
brightened and peaked before the optical /NIR synchrotron 
emission (see fig. 5 of ICoriat et all 120091 ) . suggesting that 
the jet power (which is sensitive to the high energy emission) 
takes time (days- weeks) to build up a fter the transition from 
HS S to LHS (e.g. IRussell et al.ll2010l ). 

iKalemci et"aH (|2006l ) also argued that the compact jet 
did not turn back on until well after the transition back 
to the intermediate state, and suggested that the forma- 
tion of the compact jet may req u ire th e source to be in 
the LHS. However, iFender et all (|2009l) noted the radio 
reactivation during th e 2003 outburst of XTE J 1720-318 
|Brocksopp et al.ll2005l ). which occurred marginally prior to 
the source reaching the hardest spectral state. Our data 
provide better constraints on the reverse transition, clearly 
showing a radio rebrightening between the last observation 
in the HSS and the first in the HIMS, 6-9 d before the source 
reached the canonical LHS. 

Since we detect radio emission in the HSS which we at- 
tribute to external shocks as the jets interact with the sur- 
rounding medium, we must address the question of whether 
this rebrightening is due to further external shocks, or to 
the reignition of the compact core jet. From Table [T] we see 
that the error bars on the radio spectral indices are too large 
to accurately discriminate between a flat spectrum compact 
jet and optically thin ejecta. However, since there is no fur- 
ther jump in radio brightness between the first detection in 
the HIMS and the LHS, when we believe the compact jet 
should be present, this suggests that the radio brightening 
seen between the HSS and the HIMS is indeed due to the 
reactivation of the core jet. 

The compact core jet is not detected in the individual 
VLBA data sets during the reverse transition. Stacking the 
VLBA data from the HIMS and LHS during the reverse 
transition shows a 4.5a detection (0.36 ± 0.08 mjy) close to 
the measured position of the core (see Fig. [1} . However, the 
marginal nature of the detection and the positional offset 
(22 ±3 mas in PA -14.8 ±2.9° E of N) imply that we cannot 
reliably associate this source with the compact core jet. The 
proper motion of H1743-322 is unknown, but from a com- 
parison of our measured position with that determined from 
VLBA observations made in 2004, it is unlikely to be large 
enough to have shifted the X-ray binary position 22 mas in 



a month. Since the jet axis seen by the VLBA at the peak of 
the outb urst agrees so well w ith that observed several years 
earlier bv lCorbel et al.1 (|2005l ). we can rule out precession of 
the jets. Combined with the marginal nature of the detec- 
tion, its position angle (almost perpendicular to the jet axis, 
in the absence of significant systemic proper motion) makes 
it probable that this detection is simply a noise spike. 

Given the noise level of 0.08 mjy beam -1 in our stacked 
image, we might have expected to detect the compact core 
jet, given the flux densities of 0.41-0.63 mjy measured by 
the VLA over this period. Our VLBA non-detection might 
therefore suggest that the emission measured by the VLA 
was sufficiently extended to be resolved out by the VLBA, 
and might therefore be attributed to diffuse emission down- 
stream in the jets. However, the phase decorrelation of up to 
50 per cent reported in S ect ion \3 . 2 1 would reduce the source 
flux density below the detection threshold in the stacked 
VLBA image, so the VLBA non-detection cannot be used 
as evidence for extended emission downstream in the jets. 

Assuming, therefore, that we are indeed seeing the re- 
activation of the compact core jet at the transition from 
the HSS to the HIMS, this re-ignition of the radio emis- 
sion occurs well before the optical/NIR jet is believed to 
switch on at the transition to the LHS. A simi lar sequence 
of events was also seen by ICoriat et~all (|2009l ) during the 
soft-to-hard transition in GX339-4. We can contrast this 
with the known behav i our d uring the hard-to-soft transi- 
tion. [Russeir&jiindei] (|2007l ) showed that the optical/NIR 
jets are suppressed at the top-right hand corner of the HID, 
as the X-ray emission begins to soften, while the radio flar- 
ing event and compact jet quenching occurs 7-12 d later, at 
the hypothesised 'jet line'. While we have no optical/NIR 
coverage of the 2009 event in H1743-322, the radio quench- 
ing and subsequent jet ejection did occur part-way through 
the X-ray spectral softening (see Fig. [5]l . 

Thus, it seems that the quenching and reactivation of 
the jets at the beginning and end of an outburst are grad- 
ual processes, occurring at different times for the radio and 
optical/NIR jets. Any discussion of these phenomena should 
therefore be very clear as to the wavelength at which the jets 
are being observed. A possible explanation for the observed 
sequence of events could be that the jet power evolves over 
the course of several days, with the spectral break from op- 
tically thick to thin moving to lower/higher frequencies with 
time at the beginning/end of the outburst (the frequency of 
the spectral break determines the total radiative power of 
the jets). An alternative explanation could be that it takes 
time for the suppression of the flow in the inner regions to 
propagate outwards to the surface of optical depth unity at 
radio frequencies. However, since the compact, hard-state 
jets are not resolved in VLBI observations (see Fig. [TJ , this 
would imply bulk flow speeds below 0.1c, hence we do not 
favour this explanation. 



6.2 Distance 



iMaccaronej (|2003l ) has shown that the reverse transition 
from HSS to LHS occurs at a fixed fraction (0.01-0.04) 
of the Eddingt o n lum inosity, Z/Edd- A more detailed study 
bv lDunn et all (|2010D showed that when using the fraction 
of the bolometric luminosity in the power-law component 
to identify the state transitions, the reverse transition oc- 
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curred at 0.005-0. 11/Edd , with a mean of 0.0281/Edd- As 
shown in Section [TJ the reverse transition for H1743-322 oc- 
curred at approximately the same 2.1-18.3 keV X-ray inten- 
sity in both the 2003 and 2009 outbursts. This level corre- 
sponds to a 2-20 keV fl ux of 1.5 x 10~ 9 erg cm " 2 s' 1 (ob- 
servation number 158 of iMcClintock et alj|2008h . Using the 
tabulated photon index of 2.1 appropriate for that obser- 
vation, we calculate the 0.5 keV-l OMeV flux (an a pproxi- 
mation to the bolometric flux; see M accaron 3 120031 ) to be 



5.2 x 10 9 erg cm 2 s 1 . Assuming that this corresponds to 
0.028±J^ 9 L E dd (|Dunn et alj|20ld ). we derive a distance 
to the source of 7.6l3.g(MBH/10M" Q ) 1/2 kpc, where M BH 
is the black hole mass. Using instead the range derived by 
iMaccarond (|2003l ) gives 6.5l^(M BH /10M Q ) 1/2 kpc. While 
relatively poorly constrained, the se estimates a r e cons istent 
with the distance measured by ISteiner et al.l (|2012T l. and 
thus imply a relatively high-mass compact object. 



7 COMPARISON WITH THE 2003 OUTBURST 

Fig. [7] shows the HIDs for both the 2003 and 2009 out- 
bursts of H1743-322 plotted on the same scale, with mark- 
ers corresponding to the measured radio flux densities over- 
laid (the radio data from 2003 have been taken from 
IMcClintock et al.|[2009h . The 2003 outburst was significantly 
brighter in both radio and X-ray bands, and reached a 
much softer, more disc-dominated X-ray spectrum in the 
HSS. It was a more complex outburst, with several distinct 
peaks in both the radio and X-ray light curves suggesting 
multi ple ejection events. From their model fitting of VLA 
data. ISteiner et all (|2012l ) determined the ejection date of 
the jets responsible for the downstream flaring later de- 
tected bv lCorbel et al.1 (|2005T l. which they found to be MJD 
52767. 6±1.1. Both this ejection event and the peak radio lu- 
minosity (occurring 30 days earlier, and presumably caused 
by a previous ejection) occurred at different spectral hard- 
ness values compared to the 2009 outburst, suggesting that 
the so-called 'jet line' is not constant between outbursts. The 
brightest radio emission corresponds to a harder spectrum in 
2003 than in 2009. We see no evidence for quenching of the 
radio jets prior to the flare in the 2003 outburst, although 
the cadence of the observations was not sufficiently rapid to 
rule out having missed such an event. 

While the hard to soft transition occured at a much 
higher 2-18 keV X-ray intensity in 2003, the reverse tran- 
sition occured at a very similar 2-18 keV intensity in the 
two outbursts. We also see that the radio emission as- 
sumed to arise from the compact core jet reactivated at a 
similar point in both outbursts. This reactivation occurred 
slightly prior to the hardest spect r al sta te being reached, 
in what is defined by Motta et ail (|2010h as the HIMS, or 
by IMcClintock et all ((20091 ) as a hybrid state with charac- 
teristics of both the hard and steep power-law states. As 
discussed bv lBellonil (|201Ch , such a hybrid state should cor- 
respond roughly to the HIMS. 

The HRD for the 2003 outburst is shown in Fig. gj In 
contrast to the 2009 outburst (Fig. [6} , the peak radio flux 
density occured before the sharp drop in the integrated frac- 
tional rms variability. While we do not know the exact ejec- 
tion date corresponding to the radio peak in 2003, owing to 
a lack of high angular resolution radio data, we know from 
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Figure 8. HRD for the 2003 outburst of H1743-322. The size and 
color of the markers show the radio flux density. X-ray spectral 
and timing information at the time of each radio observation has 
been interpolated from RXTE PCA ob servations within 2 d. Tim - 
ing information has been taken from IMcClintock et all l|2009t ) . 
The large sta r repre sents the moment of ejection as estimated by 
ISteiner et al with the smaller stars showing the extent of 

the la error bars on this date. Despite the location of the ejection 
event in the HRD, it occurred as the variability was increasing; 
the next significant drop in the rms did not occur until 17 days 
later. 



the measurements of ICorbel et al.l ()2005h that the jet speed 
was higher (implying a smaller delay between ejection and 
the formation of internal shocks should the core jet have a 
similar Lorentz factor in the LHS in both cases). Assuming 
that the radio peak occurred after the initial ejection event, 
the date of the first ejection in the 2003 outburst must have 
occurred before the initial drop in fractional rms. 

As noted above, the 2003 outburst was a compli- 
cated event, with m ultiple peaks in the radio light curve. 
ISteiner et all (|2012T l determined the ejection date associ- 
ated with the third major radio peak of t he outburst. When 
compa red with the results presented by IMcClintock et all 
l|2009h . their derived ejection date follows a spike in the X- 
ray count rate, and is simultaneous with an increase in the 
fractional rms variability from < 5 per cent to ~ 15 per 
cent, in stark contrast to the 2009 outburst, when we infer 
the moment of ejection to have been immediately before (or 
during) a drop in fractional rms variability (bearing in mind 
the 3-day gap in the X-ray coverage). This increase in frac- 
tional rms variability also coincide d with an increase in th e 
QPO amplitude and quality factor (IMcClintock et al"]|2009l '). 
i.e. with the clear appearan ce of classical Type C QPOs. If 
the ejection date derived bv lSteiner et al.1 (|2012T l is correct, 
then there is no unique X-ray variability signature that can 
be associated with the ejection events in both 2003 and 2009. 



8 DISCUSSION 

The 2009 outburst of H 1743-322 was a relatively simple 
event, and the lack of multiple overlapping flares provided 
a clean data set with which to test the phenomenological 
mod el for the disc-jet coupling in black hole X-ray bina- 
ries (I Fender et al.|[2009l ) . The high angular resolution VLB A 
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Figure 7. Comparison of the hardness-intensity diagrams for the 2003 (dashed line ) and 2009 (solid line) outbursts of H1743-322. 
Circles indicate radio detections (with measurements from the 2003 outburst taken from[McCliiitock ct al. 2003), with the size and shade 
of the circles indicating the measured radio flux density. Open triangles indicate non-detections. The X-ray state at the time of the 
radio observations has been interpolated from RXTE P CA observations with in 2d. The large star represents the moment of ejection as 
estimated in Section 13.3.21 for the 2009 outburst and bv lSteiner et all (|2012h for one ejection event during the 2003 outburst, with the 
smaller stars showing the extent of the la error bars on these times. The outburst in 2009 did not reach as low a hardness or as high 
an intensity as that in 2003. The ejection markers clearly show that the jet line is not at a fixed hardness. The hardness at which the 
ejection event occurs varies between outbursts. 



data allowed us to track back the radio ejecta to find the 
epoch at which the transient jets were launched, such that 
we do not need to rely on the epoch of peak radio flux 
density as a proxy for the ejection date. We now use this 
inform ation to reconsider the unified model of iFender et al.l 
attempting to answer some of the questions raised by 
IFender et ail (|2009h . However, we note the caveat that this 
is a single outburst in a single source, and that more multi- 
wavelength campaigns incorporating high-cadence, high an- 
gular resolution VLBI observations are required to verify our 
conclusions. 

Within the uncertainties created by the RXTE sam- 
pling, we find that the ejection date is consistent with the 
transition from HIMS to SIMS, at the point where the Type 
C QPOs disappear and the integrated fractional rms vari- 
ability begins to decrease. This is shortly before the detec- 
tion of quenched, optically-thin radio emission. 



The most detailed comparison to date of the connection 
between radio flaring and X- ray timing signatures was pre- 
sented by I Fender et al. | (|2009l ) (see their figs. 4 and 5). They 
investigated possible connections between the radio flaring 
events in XTE J1550-564, XTE J1859+226 and GX339-4 
and the contemporaneous drops in rms variability and ap- 
pearance of Type B QPOs. Since the exact sequence of 
events appeared to differ between the three sources, they 
concluded that there appeared to be no causal connection 
between timing signatures from the inner disc and jet ejec- 
tion events. However, they noted that any intrinsic signature 
could be blurred by both the phase of jet instability prior 
to a radio flare and the time delays incurred in producing 
optically thin radio emission from internal shocks. In light 
of the conclusions drawn from our study of the 2009 out- 
burst of H1743-322, we re-examine the flares in their three 
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sources, considering also the 2010 flare of XTE J1752-223 
and the 2003 flare of H1743-322. 



8.1 The 1998 outburst of XTE J1550-564 

The X-ray emission from the 1998 flare of XTE J 1550- 
564 peaked on MJD 51076.0 (jSobczak et all |2000| ). fol- 
lowed by a radio flare pe aking at 843 MHz on MJD 51078 
(|Hannikainen et "all l200lf) . If the radio flare followed the 
canonical pattern of an adiabatically expanding synchrotron 
event, we can infer that the peak at GHz frequencies might 
have occurred slightly prior to this. VLBI observations on 
MJD 51080.7 and 51081.7 detected several radio-emitting 
components. Assuming that the two outermost components 
in the two images can be associated with one another, the 
rate of increase in angular separation was ~ 115masd _1 , 
implying a zero-separation date of MJD 51078.6, coincident 
with the radio flare but well after the X-ray peak. The ap- 
pearance of a third component in the second image suggests 
that there could have been multiple ejection events during 
this outburst. Considering the additional uncertainties in in- 
terpreting the images (not being phase-referenced, the true 
alignment of the two VLBI images is unknown), we cannot 
unequivocally relate the timing of the X-ray flare and the 
first radio ejection even t . 

iFender et all (|200g| ) show that the radio flare began as 
the rms variability began to drop sharply. However, we also 
note a slight change in the characteristics of the Type C 
QPOs at the time of the X-ray flare (between MJD 51076.8 
an d 51077.9), moving fro m Type C to Type C' (as defined 
by iRemillard et al. 1 12002| . with weak harmonics, a double- 
peaked power density and rapidly shifting QPO frequency). 
The QPOs reverted to Type C on MJD 51078.1, and the 
first Type B QPO was not observed until four weeks later, 
on MJD 51106.95. 



steep drop in rms variability between MJD 52410.5 and 
52411.6. Again however, we remark on the change in char- 
acter of the Type C QPOs aro u nd th e time of the radio 
flare. According to lBelloni et alj l |2005l ). the Type C QPOs 
changed to Type C* between MJD 52406.7 and 52410.5, be- 
fore the first obse rvation of a Type B QPO on MJD 52411.6. 
Type C* QPOs jCasella et all 120041 ') have a broader peak, 
a higher frequency, and a lower rms amplitude than canon- 
ical Type C QPOs. Thus while the Type B QPOs and the 
drop in rms variability appear to occur after the radio flare, 
it is plausible that the evolution away from the canonical 
Type C QPOs marks the disruption of the steady jet and 
the beginning of the ejection event. 



8.4 The 2010 outburst of XTE J1752-223 

Although VLBI monitoring of XTE J1752-223 was carried 
out during its 2010 outburst, the ob served one-sided ejec- 



tion event jMiller- Jones et aL )2011bT ) and the deceleration 
of the ejecta ( Yang et alj|20ich prevent us from accurately 
determining a zero-separation time. However, the radio light 
curves and spectra show evidence for a quenching of the 
flat-spectrum jet between MJD 55217 and 55220, with the 
peak of the initial radio flare occurring on MJD 55221. This 
radio quenching corresponds to a transition from Type C 
to Type A/B QPOs betw een MJD 55217.9 and 55218.8 
(|Shaposhnikov et alJl20T(il ). In the absence of VLBI moni- 
toring close to the state transition, we cannot unequivocally 
relate the timing of the radio ejection event with the changes 
in the X-ray timing behaviour of the source. However, it is 
plausible that the sequence of events was similar to that 
which we observed in H1743-322, namely a jet ejection im- 
mediately preceding the radio quenching, and contempora- 
neous with the change in the nature of the QPOs. 



8.2 The 1999 outburst of XTE J1859+226 

The 1999 flaring event of XTE J1859+226 comprised mul- 
tiple radio flares, with the first (and largest) preceded by 
a quench on MJD 514 67.5 before peaking o n MJD 51468 
Brocksopp et alj|2002j ). The X-ray study of ICasella et al.l 



20041 ) shows a transition from Type C to Type A QPOs 



between MJD 51466.9 and 51467.6, coinciding with a steep 
drop in the rms variability (from 16.0 to 2.8 per cent) and 
also w ith the radio quenching event seen bv lBrocksopp et al.l 
|2002T ). Prior to this transition, Type C QPOs had been a 
persistent feature of the power spectrum for at least four 
days. This sequence of events appears to tie in well with 
what we see in H1743-322; a radio quenching coincident with 
the disappearance of Type C QPOs and a steep drop in the 
rms variability. If the two events are indeed similar, we can 
infer that the radio ejection event in XTE J1859+226 likely 
occurred at this time, with the short delay before the onset 
of the radio flare arising from the time taken for internal 
shocks to form and become optically thin. 



8.3 The 2002 outburst of GX339-4 



As noted by | Fender et ail |2009f ). the radio flare reported 
bv lGallo et all (|2004l ) o"ccurred on MJD 52408, prior to the 



8.5 The 2003 outburst of H1743-322 

As discussed in Section [7] the radio peak during the 2003 
outburst of H1743-322 occurred on MJD 52737, several 
days before the steep drop in fractional rms variability 
(occurring between M JD 527 47 and 52750, according to 
IMcClintock et alj|2009h . While iMcClintock et all l |200g ) do 
not classify the QPOs they detected, from the tabulated am- 
plitudes, quality factors and central frequencies, the transi- 
tion from Type C to Type B QPOs did not occur until at 
least MJD 52741. The steep radio spectrum observed at the 
peak of the outburst (a = —0.62 ± 0.01 between 4.86 and 
8.46 GHz) suggests that an ejection event was indeed respon- 
sible for t he peak radio emission. H owever, inspecting figures 
4 and 5 of lMcClintock et al] (|2009h . we see no evidence for a 
clear event in either the spectral or variability properties of 
the accretion flow that might have triggered this radio peak 
on 2003 April 8. Furthermore, as discus sed in Section [71 the 
subsequent ejection event considered bv lSteiner et al.l (|2012l ) 
seems to coincide with an increase in fractional rms variabil- 
ity and the re-appearance of canonical Type C QPOs. Thus, 
the radio emission from this event does not seem to fit the 
pattern seen in the other outbursts discussed in this Sec- 
tion, as it corresponds to different X-ray timing signatures 
than those that appeared to be associated with jet ejection 
events in other sources. However, we reiterate that the 2003 
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event in H1743-322 was much more complex than the 2009 
outburst of this sou rce, and the ejection event identified by 
ISteiner et all i|2012l ) did not correspond to the first peak in 
the X-ray and radio light curves. Therefore a full interpre- 
tation would also need to account for hysteresis effects (i.e. 
the time taken for the system to respond to the ongoing 
changes in mass accretion rate over the first few weeks of 
the outburst). 

8.6 Overall picture 

Jet ejection events have in the past been associated with 
the transition from the HIMS to SIMS, defined by a drop in 
the fra ctional rms variabi lity and the appearance of Type B 
QPOs jFender et al.l2009l ). Our radio monitoring of the 2009 
outburst of H1743-322 shows that the peak of the integrated 
radio light curve is not necessarily the best indicator of the 
date of the radio ejection event, likely owing to the time 
taken for shocks to form and become optically thin. When 
attempting to connect radio and X-ray behaviour that vary 
on timescales of days, this time delay can mask the signature 
of any causal connection. This highlights the crucial role of 
high-cadence, high-angular resolution VLBI monitoring dur- 
ing state transitions in order to pin down the exact sequence 
of events. 

In the 2009 outburst of H1743-322, we find that the jet 
ejection event occurred during or just before the sharp drop 
in the rms variability, coinciding (within the three-day un- 
certainty) with the disappearance of Type C QPOs and the 
quenching of the compact jet. The first detection of Type 
B QPOs occured several days later (albeit with the caveat 
of the sparsely-sampled transition) . This sequence of events 
is approximately consistent with what was seen in the out- 
bursts of XTE J1550-564, XTE J1859+226 and XTE J1752- 
223. However, the 2002 outburst of GX339-4 suggests that 
the drop in rms variabili ty is not causally co nnected to the 
radio flare, as noted by I Fender et al.l |2009l ). although the 
evolution away from the canonical Type C QPOs could still 
be associated with the radio flare. But even this appears to 
be ruled out by the 2003 outburst of H1743-322. We note 
that these latter two outbursts are the only ones in the sam- 
ple we considered that did not appear to be preceded by any 
significant radio quenching. Either the time sampling of the 
radio data was too sparse to detect a brief (< Id) period 
of quenching, or the physical properties (or, more specula- 
tively, the launching mechanisms) of the ejection events were 
different in these two outbursts. 

In summary, the spectral hardness at which the jets are 
ejected and the radio flare occurs appears to vary between 
outbursts (even for an individual system). While X-ray tim- 
ing studies may be better signatures of an ejection event, the 
exact causal sequence of events appears to vary betw een out- 
bursts, as originally deduced bv lFender et al.l (|2009l '). In the 
subset of cases where radio quenching is observed, the ejec- 
tion event seems to be connected to the point at which the 
Type C QPOs and band-limited noise disappear from the X- 
ray power density spectrum and the fractional variability be- 
gins to decrease. While we have shown that intensive VLBI 
observations have the capacity to pin down the moment of 
ejection, dedicated coverage of new outbursts with both high 
angular resolution radio monitoring and X-ray timing obser- 
vations is crucial to find out whether there is any common 



pattern to the sequence of events and hence any causal con- 
nection between the accretion and outflow phenomena. In 
this light, the end of the RX TE mission imp lies that prior 
to the advent of ASTROSAT (|Agrawalll2o"o6l 'l. such studies 
may only be carried out with XMM-Newton. 

The radio emission from the compact core jet appears 
to reactivate during the HIMS; it does not require the source 
to first reach the LHS, as appears to be necessary for the 
reactivation of the optical/NIR jets. We speculate that the 
suppression and reactivation of the jet near the beginning 
and end of an outburst, respectively, occur via a gradual 
decrease/increase in the total jet power, over the course of 
several days. If this is the case, travel-time arguments cou- 
pled with the unresolved nature of the core jets make this 
inconsistent with a picture in which the suppression is due 
to the compact core jet being switched off instantaneously at 
the base, with the spectral break moving to lower frequen- 
cies with time as the motion of the ejected material takes it 
past the surfaces of optical depth unity at each frequency. 
The instantaneous suppression scenario is also ruled out by 
the gradual evolution of the optical/NIR emission over the 
course of a few days. Since this emission comes from a re- 
gion < 1 light second downstream, suppression of the jet 
would instantaneously shut down the optical/NIR jet emis- 
sion, which is not observed. 

Reconciling our speculative scenario of gradual jet 
power evolu t ion w ith the internal shock scenario of 
iFender et all (|2004T l would require the gradual decrease in 
jet power during the hard to soft state transition to be 
coupled with an increase in either the jet velocity or out- 
flow flux, contradicting the naive assumption that jet power 
would be positively correlated with jet velocity and out- 
flow flux. Thus a later, separate ejection of faster-moving 
material might be more consistent with the observations 
(p ossibly due to the ejection of the corona, a s suggested 
by iRodriguez et al. I l2003l ; IVadawale et~aH I2003T ). This also 
fits in with what we know about the well-documented ra- 
dio quenching prior to m ajor radio outbursts in Cygnus 
X-3 (|Waltman et al.lll994h . Th e jet speed in that system 
is thought to be high (> 0.81c iMioduszewski et al.ll200ll ). 
whereas the rad i o que nching can last for up to two months 
|Waltman et alj I1994T ) , with the end of the quench phase 
being immediately followed by a major radio flare. Thus, if 
the same process is responsible for the radio quenching in 
Cygnus X-3 and the other X-ray binaries we have studied, 
the onset of the quench phase cannot be directly linked to 
the ejection event. A more plausible explanation might be a 
subsequent, higher-velocity ejection of material towards the 
end of the quench phase. 



9 CONCLUSIONS 

We have presented high-cadence radio and X-ray monitor- 
ing of the 2009 outburst of the black hole candidate X-ray 
binary system H1743-322, using data from RXTE, the VLA, 
VLB A and ATCA. The X-ray behaviour of the outburst was 
fairly standard, following the canonical path in the hardness- 
intensity diagram. 

Since the radio behaviour was fairly clean, likely con- 
sisting of a single bright ejection event, we were able to use 
our high angular resolution VLBI radio monitoring to track 
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the proper motion of the jet components and pin down the 
moment of jet ejection. The time of ejection appears to corre- 
spond, within uncertainties, to the beginning of the decrease 
in the fractional rms variability seen in the X-ray data, and 
to the disappearance of Type C quasi-periodic oscillations 
(i.e. to the transition from the hard intermediate state to 
the soft intermediate state). A quenching of the radio emis- 
sion occurs during or immediately after the initial ejection 
event, and there is a delay before the onset of the radio flare, 
likely due to the time taken for internal shocks to form in the 
outflow. This implies that the onset of the radio flare may 
not always be a good diagnostic of the moment of ejection, 
an assumption which has frequently been made in the past 
and may have hampered previous attempts to couple the 
radio and X-ray behaviour during black hole X-ray binary 
outbursts. 

The reactivation of the compact core jet at the end of 
the outburst appeared to occur during the transition from 
the high soft state to the hard intermediate state, and did 
not require the source to first reach the low hard state, as 
appears to be necessary for the detection of the compact 
jets in the optical/NIR. This suggests that the compact jet 
power may evolve gradually over a period of several days 
during its initial suppression and subsequent reactivation 
near the beginning and end of an outburst, respectively. 

A comparison between the outbursts of 2003 and 2009 
suggests that the jet speed in H1743-322 may be variable, 
and possibly positively correlated with the luminosity of the 
outburst. If the identification of the ejection event is correct 
in both outbursts, the X-ray spectral and timing signatures 
of the jet launching event are not constant between these two 
outbursts of the same source. However, the reverse transi- 
tion occurred at a very similar luminosity in both outbursts, 
from which we estimate a Galactic Centre distance and a 
relatively high-mass compact object for H1743-322. 

Comparison with previous outbu rsts appears to validate 
the conclusion of iFender et al.l (|2009l ) that there may be no 
single X-ray signature of the dramatic changes in the jet 
properties during a state transition. However, more intensive 
radio monitoring, at high angular resolution, is required to 
make a conclusive statement in this regard. 
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